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Abstract. Consider a particles performing simple, symmetric, non-intersecting 
random walks, starting at points 2(j — 1), 1 < j < a at time and ending at 
2(j — l) + c — b at time b + c. This can also be interpreted as a random rhombus 
tiling of an afec-hexagon, or as a random boxed planar partition confined to a 
rectangular box with side lengths a, b and c. The positions of the particles at 
all times gives a determinantal point process with a correlation kernel given 
in terms of the associated Hahn polynomials. In a suitable scaling limit we 
obtain non-intersecting Brownian motions which can be related to Dysons's 
Hermitian Brownian motion via a suitable transformation. 



We will consider a simple, symmetric random walks started at 2(j— 1), 1 < j < a, 
conditioned not to intersect in the time interval [0, b + c] , and end at c — b + 2(J — 1) 
at time b + c. Here a, b, c, c > 6, fixed positive integers. This model has several 
interpretations. One is as a uniform random rhombus tiling of an a&c-hexagon, i.e. 
a hexagon with side lengths a, b, c, a, b, c, see [3]. This translates directly to a dimer 
or perfect matching representation, see e.g. |14| . so it is a kind of two-dimensional 
statistical mechanics model. Another interpretation is as a boxed planar partition 
in a rectangular box with side lengths a, b and c, ^B]- The number of possible 
configurations, the partition function of the model, was computed by MacMahon, 
and is given by 



If we think of the random walks as the motion of particles, then at each time 
we have a certain particle configuration. By considering these particles at all times 
we get a discrete, finite point process. The purpose of this paper is to show that 
this is a determinantal point process and compute the correlation kernel in terms of 
the associated Hahn polynomials, ^U], The derivation is based on the general 
framework of and a variant of the orthogonal polynomial method. The main 
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result is theorem l3.1l bclow. The proof of that theorem also gives a proof of MacMa- 
hon's formula. A certain continuous scaling limit of this model, namely a fixed and 
b = c — > oo, converges to a model of non- intersecting Brownian motions all started 
at the origin and conditioned to end at the origin at time T. This Brownian motion 
model is a transformation of Dyson's Hermitian Brownian motion model. We will 
discuss these models in the next section and indicate how the correlation kernel 
can be computed in these models using Hermite polynomials and the orthogonal 
polynomial method. The result in this case is closely related to the work in p], see 
also [7]. In the last section we will consider the discrete model where the orthogonal 
polynomial method is less obvious. At the end of that section we will give some 
remarks concerning asymptotics. 

2. General framework and Dyson's Brownian motion 

2.1. General framework. Let X r , < r < to be subsets of R, (f> r , r +i ■ X r — > 
X r +i, < r < m, given functions and /i r a measure on X ri 1 < r < to, e.g. 
Lebesgue or counting measure. An element x = (x 1 , . . . , x m ~ l ) £ X™ x X% x • • • x 
X^-i = X is called a configuration. We think of x\, . . . x r n , x r = {x\, . . . , x r n ), 
as the positions of particles in X r , which we will call line r. Let x° £ Xq and 
x m 6 be fixed configurations, the initial and final configurations respectively. 
Define <fi r s : X r x X s — > R for r < s by 

(2.1) 

4> r>s (x,y) = J (f> rtr+1 (x, Zi) ... (f> a - lia (Zr- a -l,y)dfJL r+ i(zi) ■ ■ ■ d/i s _i(z r _ s _i), 

and <j) r>s = if r > s. We will consider probability measures on X of the form 

rn~ 1 

(2.2) — [] det^r+^xlx^W^dtfix 1 ) . ..dfC-iix™- 1 ), 

' r=0 

where Z n<m is a normalization constant. It is proved in |11| that the measure 
(12.21) has determinantal correlation functions, i.e. the probability density with 
respect to the reference measure d/i ri (yi) . . .dfJ,r k (Uk) of finding particles at z\ — 
(ri,yi), ...,z k = (r k ,yk) is given by 

(2.3) det(K n 

where K is the so called correlation kernel. This kernel is given by 

n 

(2.4) K n>m (r,x;s,y) = -(t) r>s (x,y) + c/) rym (x,x' j n )(A^ 1 ) l j<j) QtS (x Q j ,y), 

where A — (</>o,m(^°, ar J*))?'7=i- Note that the kernel is not unique. We can multiply 
it by ip(r,x)/i/)(s,y) for an arbitrary function ip ^ and get the same correlation 
functions. 

2.2. Dyson's Hermitian Brownian motion. Let H(t) be an n x n Hermitian 
matrix whose elements evolve according to indepenent Ornstein-Uhlenbeck pro- 
cesses, see 0], 221- We consider the stationary case. The probability measure for 
seeing the matrices Hi,. . . , -ff m _i at times ti < • •• < t m is 

1 „-tr h? TT „™ ( tr ( H i+i ~ 9jHj) 2 ■ 




(2.5) ^— e~«^ | | exp ( - " J \_£ J ' ) dH x . . . dH n 
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where dHj is the Lebesgue measure on the space of Hcrmitian matrices, and qj — 
exp(— (tj+i — tj)), 1 < j < m — 2. Integrating out the angular variables using the 
HarishChandra/Itzykson-Zuber formula, |15j . gives the eigenvalue measure 
(2.6) 

m - 2 ' ' (X^-qrXVf 



\ m — 1 \ 

where A n (A) = rii<i<j<n(^i — ^j) ^ s the Vandermonde determinant, and \ r -, 1 < 
j < n, are the eigenvalues of H r . 

If we set <j)Q y i(i,x) = pi(x)e~~ x , TO _i jm (x,i) = Pi{x), where is a polynomial 
of degree i, (j) r , r+ x(x,y) = (vr(l - qf)Y 1/2 cxp(-(y - g r x) 2 /(l - q 2 )), X = X m = 
{0, . . . ,n — 1}, irj = =z — l, X r = R, 1 < r < m and fi r the Lebesgue mea- 
sure, we see that <|2.6ll is of the form 12. 2|) . This is a basic example of a measure 
of the form 12.21) . Here we have used the classical trick in the orthogonal polyno- 
mial method in random matrix theory to write the Vandermonde determinant as 
A„(A) = det(pj(Aj)). The polynomials p, can be arbitrary but we choose them to 
be the normalized Hermite polynomials. This will lead to a formula for the kernel 
(|2.4(1 in terms of the Hermite polynomials. The key is the expansion, see e.g. pQ, 

1 ( - ) 2 00 

(2.7) e ~^r- = J2 Pk (x)p k (y)q k e-y 2 , 

< q < 1. Repeated use of this identity gives 

(2.8) Mj,y) = e-^-^ Pj (y)e-y 2 . 
Similarly, 

(2.9) 4>r, m (x,j) = e- j ^--^ Pj (x). 

Using the orthonormality we obtain (j>o,m(h j) — ex p( — j{tm-i ~ ti))Sij an d hence 
(A~ r )ij = exp(j(i m _i — ti))Sij. It also follows from (|2.7|l that if 1 < r < s < m, 
then 

I / t e U-ts x _ y Y< 

* rAX > V) = v /4l-e^)) CXP l"l-e^) 

OO 

(2.10) = Y^Pk^Pkiy^-^e-y 2 . 

k=0 

Set Xt,s = 1 if i < s and Xt,s = if t > s. From l|2.4(l we get the extended Hermite 
kernel, 

1 / (e*~ s x-y) 2 \ 

^• H ™'(^ ;S '^ = ~ v / 4l _ e 2( t -,) ) eXP I,' l_ e 2( t -s) J Xt,„ 

(2.11) +Y / e k{t ~ s) Pk(x)p k (y)e-y 2 

Using the second equality in l|2.10|l we obtain the alternative formula 

'EVoe k{t - s) Pk(x) Pk (y)e-y 2 ,t>s 
r ET= n e Ht - s) Pk(x)p k (y) e -y 2 ,t<s. 

Multiplying with exp(— x 2 /2 + y 2 /2) we get the ordinary Hermite kernel when t = s. 



(2.12) K cxt 

.Hcrm. 

(t,x;s,y) 
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Let 7 be a positively oriented circle around the origin with radius r > 0, and T 
the line R 3^ L + it with L > r. Using the integral formulas, see pQ, 

HJx) = ^-=e x2 [ e w2 - 2xw w n dw, 
W n Jr 

w 2m J 1 z n + l 
and p n (x) — (v / 7r2"n!)~ 1 / 2 iJ n (a;), it is not difficult to show that 

K eKt . HeIIU Xt,x;s,y) = -_^__exp \ - 1 _^ t _. ) j Xt,. 



(2.13) + j- .„ / tin; / e 

This double contour integral can be useful for asymptotic computations, for example 
to show convergence to the extended Airy kernel when we have the edge scaling. 
To our knowledge the details for this has not been presented in the litterature, 
but using (|2.13|) and the integral formula for the extended Airy kernel it should be 
possible to do this similarly to what was done for the extended Airy kernel in |12j . 

2.3. Non-intersecting Brownian motions. A second closely related example is 
the following which involves non-intersecting Brownian motions. Consider n non- 
intersecting Brownian motions started at x® = e(i — 1), 1 < i < n, at time and 
conditioned to end at the same points at time T. Let x\, 1 < i < n, denote the 
positions at time r r , 1 < r < to, where = tq < t\ < ■ ■ ■ < r TO _i < r m = T. 
By the Karlin-McGregor theorem the probability density for x = (x 1 , . . . , x" 1 ^ 1 ) E 
(R n ) m_1 is given by 

m—X 

(2-14) — J] det(p T ^ 1 - Tr (x r i ,x r + 1 ))Z j=1 . 

In the limit e — * 0+, corresponding to all particles starting at the origin at time 
and ending at the origin at time T, we get the probability density 

1 n — 1 m— 2 

-^(y 1 ) J] e ~W** [J det(e-K +1 -vr)V 2 (- +1 -r r)) „. =i 

j=0 r—1 
n-l 

(2.15) x A„(z/ n_1 ) J] e-ivr 1 ) 2 /^-^). 

j=o 

This has again the general form l|2.2|l with 4>o,i(i, y) — Qi{y) ex p( — 2/ 2 /2ti), (f> r ,r+i { x i y) 
exp(-(j/-.x) 2 /2(r r+ i-r,.)) and (j> m -i, m (x,i) = g. t {x ) exp( -x 2 /2(r-r m _i)), where 
qi and cji are polynomials of degree i. The measure ()2.15f) is actually a transforma- 
tion of the measure l|2.6(l . Define 



(2.16) d r = 



T 



2r r (T — r r ) ' 



1 < r < to, and r r = T(l + e 2 * r ) 1 . If we set = yjd r , 1 < j < n, 1 < r < m, 
then a straightforward computation shows that 1)2. 15[) transforms into (|2.6|) . In this 
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way we can also transform the extended Hermite kernel H2.11J) into a correlation 
kernel for Ij2.15|l . However, let us indicate how we can obtain it directly. 
Set 



/2 f r r (T-r r+1 ^ j 1 

T r +l(T — T r ) 



Then 



(2.17) / e 2 (^+i--r) e 2r r pj (^ x d r )dx = c r jpj(yd r+ i)e 2T r+i , 
Jr 

where pj is the j:th normalized Hermite polynomial. This can be deduced from the 
identity 

(2.18) f e-^ 2 p n (ax)dx = tt^I - a 2 ) n ' 2 p n { a \ ' ), 
Jr [l -a 2 ) 1 ' 2 

which in turn follows easily from the generating function for the Hermite polyno- 
mials. Choose qj(x) = pj(xd\) and q"j(x) = p 3 (xd m -i). It follows from (|2.17(l that 



1/2 

(2.19) <t> , s (j,x) = ( y l[(r i+1 - n ) ) 2(*- 1 )/ 2 Y[c hlPj (xd s )e- x2 /^ 



and 
(2.20) 

( T m-2 \ 1 / 2 m-2 

i—r / i—r 

Using the orthogonality of the pj:s and the general formula (|2.4II we obtain the 
following expression for the correlation kernel 

1 (*-y) 2 

K B M(Tr,x;T s ,y) = : e 2 (Tr - Ts) Xx)Tj 

A /27r(r r - r s ) 

(2.21) 

() K T S (T — T r ) J \2t s (T-t, 
Here we have multiplied by the unimportant factor 

\[T r ~ T r -!) . . . (T 2 - Tl) J 

If we go back to the transformation discussed above we see that 
(2.22) 

1 „ / X V s ( T I T r - l\ 1/4 

K B M{T r , -j-;t s , — )e r t — = A oxt .HGrm.(ir,a;;i s ,2/), 



\fdr~d~s dr ' d s \T/t s — 1 

with r r = T(l + cxp(— 2£ r ))~ 1 and -Kcxt.Herm. given by (|2.11|) 
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3. The extended Hahn kernel 

3.1. Derivation of the kernel. Consider a symmetric, simple random walks with 
initial points (0, 2j) and final points (b + c, c—b + 2j), < j < a — 1, conditoned not 
to intersect in the whole time interval [0, 6 + c]. The single step transition kernel 
for one particle is 

(3-1) ^(f>{x, y) = ^S x -i, y + ^S x+1 . v . 

The configuration at time t = r, which we also call the configuration on the r:th line, 
is given by points zj, < j < a, zg < ■ ■ ■ < z 7 a _ ± , where z° = 2j, z b+c = c — b + 2j. 
We think of these points as the positions of particles. By the Lindstrom-Gessel- 
Viennot method, our probability measure on the set of configurations z = (zj) 
in (Z a ) b+C ^ 1 is 

(3-2) POO = z^b-cj n^W'^ 1 ))^- 

Here Z(a, b, c) is the total number of configurations and is given by MacMahon's 
formula 1)1.1)1 , 

The measure (|3.2ll has exactly the general form (|2.2|l (with /x counting measure 
on Z), and we want to compute the correlation kernel 1)2. 4[l . To do this we will 
use the orthogonal polynomial method in a similar way that was used for the non- 
intersecting Brownian motions in the last section. How this should be done is not 
obvious from (|3.2|l . It is shown in [5], that the induced probability ensemble on a 
single line is an orthogonal polynomial ensemble, where the relevant polynomials 
are the associated Hahn polynomials. This indicates that we should modify the 
first and the last factors in 1)3.2)1 by doing row operations so that we get a situation 
where the matrix A in 1)2.4)1 is diagonal. 

The normalized associated Hahn polynomials, \Ijg\, HO]- 0- can be defined using 
a hypergeometric function by 

(a.ffi, s (-JV -/?)„(- JV) n {- n ,n-2N -a- 0-1,-x , 

q n,N^)- -J^T- 3 2 I -N-p-N ^ 

a n,N V 

(-JV - p) n (-N)n A fn\ , v , {-x) 3 { n -2N-a-/3-l)j 



(3 - 3) ~ d^fnl l) {-N-PU-N), 

where 

/„ 4) (A^0)\ 2 = (a + (3 + N + l-n) N+1 

1 ' ' \ n > N J (a + (3 + 2N+l- 2n)n!(/3 + JV - n)!(a + N — n)!(JV - n)\ ' 

and we use the standard notation (a) n = a(a+l) . .. (a + n— 1). These polynomials 
are orthogonal with respect to the weight 

(3 ' 5) WN^ix) = x \( x + a )\( N + - X )\( N - x )\> 

on {0,1,..., JV}, i.e. 

N 

(3-6) T,&\^&\^> { N' P) ^) = 5n,m, 

x=0 
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for < n, m < N. Below we will sometimes use the convention that 1/nl = if 
n < 0, so that the summation in Ij3.6|l for example could be extended to x G Z. 

Our goal is to give a formula for the correlation kernel in terms of the associated 
Hahn polynomials. First, we need some notation. Let a, b, c G Z + , b < c. Set 
a.,. = \c-r\, b r = \b-r\, 



(3.7) 
and 



-r , < r < b 
r -2b ,b<r <b + c 



{r + a-1 ,0 < r < 6 

b + a-1 ,b<r <c 

a + b + c — 1 — r , c < r < b + c. 

Define 

{{(b r + x)\(j r + a r - x)!)- 1 ,0<r<6 
(x!(7 r +a r -x)!)" 1 ,b<r<c 
(x\(j r - xy.y 1 ,c<r<b + c 

and 

Uylhs-yVr 1 ,0<s<b 
(3.10) u s (x) = I ((6 s + y )!( 7s - 2/ )!)-i ,b<s<c 

{ ((b s + j/)!( 7a + a s - y)!)" 1 ,c<s<b + c. 

Theorem 3.1. The point process on {1 a ) b+c ~ 1 defined by \3.21\l has determinantal 
correlation functions with kernel given by 

K H (r,a r + 2x\ s, a s + 2y) = -0 r , s (a r + 2x, a s + 2y) 
(3.11) 



a-l 



v (a + r — 1 — n)!(a + o + c — 1 — n)! 

n— V v y v y 

/or < r, s < 6 + c, ijeZ. ffere </> ns = if r > s and 
(3.12) <Mz,y) 



3/— g-j-g— r 
2 



z/r < s. 
Proof. Set 

for < j, j < a, 

and 



1 

Cj ' k ~ RmiFh)!' 



fn.k 



f* = 



n\ (n — 2a — 6 — c + 1) 



fcy (-a - c+ l)fc(-a)fc 
n\ (n — 2a — o — c + !)& 



k) (-a- b+ l)fe(-a)fc 
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for < k < n. Define 

n a — 1 

^{n,z) = ^2 / n , m ^2 c j>m <j)(2j,z), 

m— j=m 
n a — 1 

(3.13) ip*{n,z) = ^2 fn,m E Cj, m , 0(c - 6 + +2j» 

< n < a, z G Z. 

We will now do row operations to modify the first and the last factor in (|3.2|) . 
det(0(x?,4))^i o = det(0(2m,4))^ =o 

a-l 



det (r^- E c ^mm,4)T^L 

a — 1 ^ a— 1 

II ; det(^c J> <^(2i,4))^ 1 =0 



m— ' 



a — 1 ^ .j n a — 1 

n ; — det (-r- E E ^»^(2j ; 4))r fe Lo 

m— u ' m— U j—m 



a-l 



(3-14) = J] - \ det(^(n,4))^ fc L . 

Cm .mjm.m 



m=0 

In the same way we obtain 

(3.15) det(^ +c - 1 ,4+ c ))-i = H - * detWn, x^ 1 ))^. 

m=0 m ' m J m,m 

If we now set o ,i( n >2/) = tp(n,y), cj) b+c -t <b+c {y,n) = ip*{n,y) and cf> rtr+1 {x,y) 
4>(x, y), l<r<6+c— 1, the probability measure 13. 21 can be written 

(3 - i6) = WW n c2 r f= n^c^icyj.gr- 1 ))^, 

where y° = j/j +c = j, < j < a. 



Write 0*™ (a;, y) = * • • • * </>(:r, y) (n factors) if n > 2, 0* 1 (a;, y) = </>(a;, y) and 
v 



6*°(a;, y) = S x ^ y . We want to compute (j>o, s , <j> r ,b+c and </>o,6+e for 1 < r, s < b + c 



By definition 

(3.17) 0o,r(n,y) = ^^(n,z)f ^^(z.y), 
and 

(3.18) ^ r>6+c (y,r 1 )=^^(ri^)r( 6 + c - r - 1 )(0 ! y), 

zGZ 

since = (j>(y,x). 

Claim 3.2. If z e 2Z + 1, then 

(n — 2a — b — c + l)j 
,._„ w (-a-c+l)^-^^-^-^)!' 



(3.19) ^^) = E(")t^ 
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and if z £ 2Z, then ip{n 1 z) = 0. 
Proof. By definition 



(n \ (n - 2a - b - c + 1) 3 ^ 



(-a - c - l) m (-a) m ^ (a - m)(j - m)!(a - 1 - j)!' 



ri \ / V — lift \ / i it> 

m— U j—m 

Now, with 2 = 2C - 1, 

a-l 



E 



^2j-l,z + <>2i+l,z 



1 



(a - m)(C - to)! (a - 1 - C)! (a - m)(( - 1 - m)!(a - C)! 
1 



(«±i -m)!( a - *±i)!' 



2 2 

and follows. □ 

If z € 2Z — 1, then simailarly 

(rt — 2a — o — c + l)j 



(3.20) <r(n,c-6 + *)=X;Q 



( _ a _6 + i).(_ a ). ( *±i_j) !(a _»±i)r 



Claim 3.3. 

(3.21) 



. f n\ in - 2a — b — c + 1) , 

, r (n,y) = E Jt-T— XTT737— X 



£5 U7 (-a - c+ lj^-a - r + 1), _ j)!(a _ i _ 2/_!)i " 
Proof. By induction on r. The statement is true for r = 1 by l|3.19|l . We have 

^ o ,r("-,a;)(4, a +i + 5x,y-i) = 4>o,r(n, y + 1) + <fo, r (n, y - 1) 

'ra\ (n- 2a- b- c+ l)j 1 



x a 



^ W (-« - c + l)j(-o - r + 1), (I±I±1 - j)!( a - 1 - ^)! 
y — r + 1 y + r + 1 



2 2 
Now, 



j1- 



(o + l) r -l , ■ \ (a+l)r 

-(a-r-j) - 



(-a - r + l)j ' (-a - 

and the claim is proved. □ 

Also, 

0r,6+c(c - 6 + x, n) = J2 ^*("« 2)r (b+c ~ r ~ 1} (^ c - 6 + x) 

zGZ 

= ^ ^*(n, c - 6 + z))0* (b+c -'- 1) (c - 6 + z, c - b + x) 

== E ^*^ n ' c - b + 2))0* (b+c_r_l) (z, x). 
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We can now proceed exactly as in the proof of claim EP1 and show that 
, \ , ,s v^AA (n — 2a—b — c+l)j 

(/)r.b+c{y,n) = (a+ 1 b+c-r-i > . 7— — ; ; — r- 

(3.22) 1 



(^ + b-j)l(a + c-l-y±L)\- 

Introduce new coordinates, which we will call the Hahn coordinates on line r by 

_ y r k -ar 
Xk ~ 2 ' 

Then, < x r k < j r . One motvation to use these coordinates is that it is easier to 
recognize the Hahn polynomials when using them. Since (j>o tr (i,z) is zero unless 
z + r is even, i.e. unless z — a r is even, we obtain 

(3.23) A nm = 2J 4>o,r(n, a r + 2z)4> ryb+c (a r + 2z, to). 

The correlation kernel is given by 

K{r, 2x + a r \ s, 2y + a s ) = -4> r , s (2x + a r , 2y + a s ) 
a-l 

(3.24) + (l>r,b+c(2x + a r ,i){A- 1 

)ij ( t ) o,s(j, 2y + a s ) 

i,j=0 

according to l|2.4[) . We want to express (f>o jr (j,2y + a r ) and 4> r fi+ c {2x + a r ,i) in 
terms of the associated Hahn polynomials. In order to do so we have to distinguish 
three cases, 1 < r < b, b < r < c and c < r < b + c. 

Set a r = |c — r\ and b r = \b — r\. 
(i) Consider first the case 1 < r < b. By $F3$ and (|?HT|) 

. AA (n — 2a—b — c+l)j 
^o, r (n, a r + 2z = a + 1 r -i V . y-^ -777 7-77- 

(3.25) 

v 1 _ (a + l) r -idl b ; 7 ; a '' ) n! (6r<ar) 



(z-j)!(o + r-l-«)! (-a-c+l) n (-a-r + l) n ^ w z!( 7r -z)f 
Also, by E22J), 

(n-2a-6-c+l) 3 - 



j7 (-a-i»+l)j(-a-&-c + r + l)j 



r,HcK + 2z i n) = (a + l)b +c -r-l 

3=0 
1 

(6 - r + z - j)\(a + c-l- z)\ 

(a + l) b+c _ r _ 1 ^ (—n)j(n — 2a—b — c)j(—b + r — z)j 



E 



(b-r + z)\(a + c - 1 - z)\ j\(-a-b+ l)j(-a -b-c + r + l)j 
(a + l)b+c-r-i p {—n.n—2a — b — c + l,—b + r — z ^ 



(6,. + z)!(7 r + a r - z)\ V -o-6+l,-a-6-c + r + l 

We can rewrite this using the following hypergeometric identity, yQ p. 141, 
(3.26) 

f-n,a,b \ _ (d—a) n (e — a) n f-n, a, a+ b - n - d - e + 1 x 
3 2 { d,e ; J- (d)„(e) n 3 2 ^a-n-d+l,a-n-e+l 5 / 
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This gives 

(a + l) 6+c _ r _ 1 (a + c - n)„(a + r - n) r 



<f>r,b+c(ot r + 2z) 



(-a + b+ l)„(-a - 6 - c + r + l)„(b r + z)!(7,, + a r - z)! 

/-n, n - 2-/ r - a r - b r - 1, -z 
x 3^2 ; i 

(a + l)fc +c _ T ._i(a + c - rt) n (a + r - n) n dn, r ^n\ 

(-a + 6 + l) n (-a - b- c + r + 1)„(— a - c + l)„(-a - r + 1)„ 

(3.27) 

X ^ ;ar)W (6, + z)!( 7r + ar -z)! 

We can now compute A nm given by (|3.23|) by picking r between 1 and b, the choice 
does not matter. Using i|3.6[l . 1|3.25[1 and i|3.27[l we obtain, after some simplification 

(3.28) A nm = C n (a, b, c)~~ 1 5 n , m , 
where 

(a + 6- l)l(o + c - l)!(2a + H-c-2n - l)a! 2 

(3.29) C n (a, 6, c) = . 

(ii) Next we consider the case b < r < c. The computations are similar to those in 
the previous case. We find 



4>Q,r{n, a r + 2z) 



(a + l) r _i(a + b — n) n (a + b + c — r — ^)n^rv^" r ^n! 



(3.30) x &;^(z)- 



(—a — c + 1)„(— a — r + 1)„(— a - 6 + 1)„(— a — b - c + 1 + r) r 
1 

X (BST'W (b r + Z)l(rf r - Z)V 

Here we have used the hypergeometric identity i|3.26[) . Also, we find 

(a + l)b+c-r-\dn^ T ' '«! 



^f, +c (a r + 2z,n) = 



(-a - 6 + 1)„(— a - 6 - c + 1 + r)„ 
1 



(3-31) xq^(z) — 

z\(j r + a r — z)\ 

(iii) Finally we come to the case c < r < b + c, and again the computations are 
similar. We obtain 

(a + l) r _i(a + b — n) n (a + b + c — r — n) n dn^ r ^n\ 



</>o, r («, a r + 2z) = 



(-a - c + l)„(-a - r + l)n(— a - 6 + l)n(— a - 6 - c + 1 + r)„ 



(3 - 32) X ^^^^ (6 r + z )\{ lr + a r - z)V 

where we have used the identity l|3.26|) . Also, 



4> r ^ +c (a r + 2z,n) 



(a + l)b+c-r-idn^ r ' ] n\ 
(-a-b+ l) n (-a - b - c+ 1 + r) r 



(3-33) x ^) (a) _J_. 

We now have all the ingredients in (|2.4p . It follows from (|3.28|l that 

(A" 1 )^ = Ci(a,b,c)5ij, 
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and some computation now gives (|3.11|l . Note that 4> r , s {x,y) is the number of 
random walk paths from x to y in s — r steps and hence is given by Ij3.12|l . 

The computations in the proof of the theorem also gives a proof of MacMahons 
formula. We have 

a— 1 ^ 

(3.34) Z(o,6,c) = n 2 a a- detA 



n=0 



A computation gives 

jj 1 _ jj (2a + b + c-2n- l)! 2 (a + b~ l)!(a + c - 1)! 

11 r2 w„ „w* 11 



cl,Jn,nd* n>n ^ (a + & - 1 - n)!(a + c - 1 - n)!(2a + 6 + c - n - l)! 2 ' 
It follows from (|3~3g|) and (|3~^|) that 

TT n\(2a + b + c-n-l)\ 

11 (a + 6- l)!(a + c- l)!(2a + fc + c-2n- l)a! 2 ' 

Hence, by l|3.34|l and after some simplification 

Z(a,b,c) =[[ 



■J- (b + n)\(c + n)V 



n=0 



which is the same as (|l.lf) . □ 



3.2. Some remarks about asymptotics. As discussed above the non-intersecting 
Brownian motion model 12.15fl is a kind of continuum version of the random walk 
model. In fact it can be obtained as a scaling limit of the random walk model. For 
the associated Hahn polynomials we have the asymptotics 
(3.35) 

lim d£&n\ | ; 2 = I p£g>(- + 2zJ----N) = H n (z), 

N^oc n < N y Af3/2^(2t + l)(i + l) J p ^ Ny 2 V t + 1 ' nW ' 

where a/N — > t > 0, uniformly for z in a compact subset of C. Here H n (z) 
is the ordinary Hermite polynomial of degree n. This can be proved by a slight 
modification of the argument in |S] based on the recurrence relation. Using (|3.35|) 
and standard asymptotics for the binomial coefficient it follows that 

(3.36) 2 r - s ^K H (r, x; s, y) - K BU (r, £ a, y) 

as A: — > oo if r/fc -> 2r/T, s/jfe -> 2ct/T, ir/Vfc -> S.y/2/T, y/Vk -► rj- s /2/T, where 
Kbm is given by l|2.21|l . So in this sense we have convergence to the Brownian 
motion model. It should also be possible to prove this directly, i.e. that the 
measure <|3.2[1 converges, when rescaled as above, to the measure (|2.15|) . compare 
the arguments in |13| . 

A more interesting, and also much more difficult limit is to consider the case 
when, a, b and c go to infinity with the same rate, say a = b = c — > oo. In 
particular it is interesting to consider the fluctuations of the top (and bottom) 
curves which bound the so called frozen regions, [3], in the tiling. If we restrict 
to a single line, this has been done recently by |2] using very precise asymptotics 
for Hahn polynomials derived using Ricmann-Hilbert techniques. This shows for 
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example that if n = a = b = c then the last (first) particle fluctuates like n 1 / 3 
in the appropriate region and that the fluctuations are given by the Tracy- Widom 
distribution, see [2] for details. If these asymptotic results could be extended to 
the extended (associated) Hahn kernel, l|3.11|l . it should be possible to prove the 
convergence of the boundary curve of the frozen region to the Airy process, |17| . 
|1 1 j - as has been done for some other tiling problems in and |12j . 
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